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►^ , Abstract 

O 
^ . In a previous paper an asymptotic expansion for A^ in powers of 1/d was developed. 

The results of computer computations for some terms in the expansion, as well as vari- 
'nI" ■ 

O ' ous quantities associated to the expansion, are herein presented. The computations (in 

OO . 
O ■ integer arithmetic) are actually done only for d = l,d = 2, and d = 3, but some results 

for arbitrary dimension follow from the general structure of the theory. In particular we 

rS I obtain the next term in powers of 1/d, the l/rf^ term in the asymptotic expansion for A^, 

as well as checking the correctness of the terms previously calculated by hand. 



For d = 2 and d = 3 a. second expansion for A^ is defined and numerically studied. 
Contrary to our hopes this second type of expansion appears not to be convergent. There 
is not enough evidence to argue if the expansions for A^ in powers of 1/d are convergent 
for either d = 2 or d = 3, but likely they are not. ( Certainly convergence, but at a very 
slow rate, remains an interesting possibility in all our cases. ) 



In [1] an expansion for Xd in powers of 1/d was presented, completing the asymptotic 
expression obtained in [2], by Mine. The expression as presented in [1] is 

A.~-M2d)-- + -- + --p+... (1) 

One resuh of this paper is the extension for one more power of 1/d 

1 , / n 1 1 1 5 1 5 1 
A.~-te(2<i)-- + -- + --^ + --^ + .... (2) 

The key quantities to be computed are the Jj defined in equation (28) of [1], and the 
note after (33) therein 

j.= j^ (3) 

The first few Jj are given in (29) and (38) of [1]. 

We have used Maple software to compute the Ji for i = 1,2,3,4,5,6; d = 1,2,3 as 
given in the following. 



(4) 
(5) 
(6) 
(7) 
(8) 
(9) 

(10) 

(11) 
(12) 

(13) 

(14) 

(15) 



For d = 1 one has: 










Ji 


= 







J2 


= 


1/8 




J3 


= 


1/12 




Ja 


= 


-3/64 




J, 


= 


- 13/80 




Je 


= 


- 19/192 


For d = 2 one has: 










Ji 


= 







J2 


= 


1/16 




J3 


= 


1/48 




J4 


= 


- 9/512 




h 


= 


- 23/1280 




Je 


= 


25/3072 



For (j = 3 one has : 

Ji = (16) 

J2 = 1/24 (17) 

J3 = 1/108 (18) 

J4 = - 5/576 (19) 

J5 = - 11/2160 (20) 

Je = 175/46656 (21) 

The resuhs above for Jj, i < 4, were presented in [1] as obtained by hand computation. 

The Ji for d = 1 are computable in seconds of computer time, and Jj with much 
higher values of i are easily obtained. For d = 2 and d = 3 the Jq required several days 
of computer time; J7 would be extremely difficult to obtain, requiring both more theory 
and much more computer time, a year's effort. 

The results, (4) through (21), and equation (30) of [1] (whose proof will appear in a 
later paper), will establish the following formulae for Ji in general dimension: 

Ji = (22) 

7 _ 3 1 3 1 

'^^ - ~32rf^+64rf^ ^^^^ 

J = ll_ll_il ^271 

^ 48ci3 64^4 5^5 y ^ 

Substituting these values of Ji into (36) of [1], and extracting the behavior of the 
exponent in (36) to order 1/d^, we get the result in (2) above. We are using the fact that 
Ji for i > 6 have only powers of 1/d greater than 3. 



For d = 2 and c? = 3 we have considered two sequences of approximations each 
characterized by keeping more and more terms in some sense. The first sequence of 
approximat ions 

Ao,^i,--- (28) 



A.^iM2.)-i + ^ + | + ... (29) 



Ar = l ln{2d) - i + ^ |. (30) 

1=1 



A, = 


-- .1931 


Ai = 


= .2556 


A2 = 


= .2687 


As = 


= .2784 



is constructed from 
by 

We then have the results: 

For d = 2 : 

(31) 
(32) 
(33) 
(34) 

Of course here, for d = 2, the exact result A2 = .29156.... is known. [3], [4] 

For d = 3 : 

(35) 
(36) 
(37) 
(38) 

Exact bounds for A3 as given in [5] are 

.440075 < A3 < .457547 (39) 

Our second sequence of approximations is 

5o,fii,--- (40) 



Ao = 


-- .3959 


Ai = 


= .4375 


A2 = 


= .4433 


As = 


= .4462 



constructed by replacing Ji by Jj a;*"^ with x a formal parameter. We then substitute 
these Ji x^~^ for Jj in (36) of [1]. To obtain B^., keeping terms through power x^ (and 
then setting x = 1), we get the values: 



for d = 2: 



Bo = .1931 (41) 

Bi = .2556 (42) 

B2 = .2921 (43) 

Bs = .2992 (44) 

54 = .2905 (45) 

Br, = .2814 (46) 



and for d = 3: 



Bo = .3959 (47) 

Bi = .4375 (48) 

B2 = .4538 (49) 

B3 = .4524 (50) 

54 = .4468 (51) 

^5 = .4445 (52) 

Originally we hoped the sequence of Bi converged, but this now seems unlikely. If 
the sequences of Bi for d = 2 and d = 3 behaved similarly, we might have an argument 
that A3 = .453 ± .001. We do not know whether to believe this or not. 

It remains to outline the general flow of the computer computations. The quantities 
of interest were obtained through a chain of operations successively performed. Were 
we to redo the hand calculation of J4 as described in [1], we would flnd it much more 
efficient to echo some of the ideas used in the computer programs. Naturally we cannot 
describe the many little clever ideas, saving memory and increasing efficiency, needed to 



push the computations to the level of obtaining Jg. We find it convenient to describe our 
computations as divided into four steps, which we present in inverse order. 

Step 4. This step assumes knowledge of the Jj, and substitutes them into (36) of [1] to 
obtain the approximations Ai and -Bj. 

Step 3. This step assumes known the quantities 

XI f{si)---f{ss)'ip'c{si,S2,...,Ss) (53) 

Sl,S2,...,Ss 

(where / is as in [1], ^ on a dimer, zero on other tiles, and (53) is as (28) of [1] with v 
replaced by /) and computes 

Js= X v{si)---v{ss)i''c{si,S2,...,Ss). (54) 

Sl,S2,...,Ss 

Step 2. Takes the output of Step 1 described next to yield the expression in equation 
(53) above. 

Note that in equations (53) and (54) above the sums are over all located tiles si, ...,Ss 
such that the sequence is connected. That is, one cannot pick two subsets of ordered sets 
Si, ...,Ss such that, with the two sub-ordered sequences 

Sai ,•••, Sar and Sp^ ,..., 3/3, 

one has 

1) the same Si appear in si, ..., Sg and in the union of the two sub-ordered sets, with 
the same multiplicity 

2) the two sub-ordered sets of located tiles have no overlap between them. 

We now define an s- dimer tree to be an ordered sequence of s located tiles Si, S2, ..., s^ 
such that for each i < s the sequence of located tiles si,S2, ■■■,Si is connected. It is easy 
to see s- dimer trees are the sort of thing a computer can easily generate. 



Step 1. The output of Step 1 (towards computation of Jg) is the set of all possible s-dimer 
trees with given Si. With each s-dimer tree the output of Step 1 also includes knowledge 
of the overlap pattern of the tree. That is, for the tree Si, ..., s^ one knows whether Sj n s^ 
is empty or not for each i and j. (The output of Step 1 actually yields the number of 
s-dimer trees with each overlap pattern.) 

Steps 3 and 4 take seconds to run in the cases we have studied. In the computation 
of Jg) Step 1 and Step 2 required several weeks of computer time. 
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